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We derive the evolution equations for a system of neutrinos interacting among themselves and with 
a matter background, based upon the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. 
This theoretical framework gives an (unclosed) set of first-order coupled integro-differential equa- 
tions governing the evolution of the reduced density matrices. By employing the hierarchy, we first 
rederive the mean-field evolution equations for the neutrino one-body density matrix associated with 
a system of neutrinos and anti-neutrinos interacting with matter and with an anisotropic neutrino 
background. Then, we derive extended evolution equations to determine neutrino flavour conversion 
beyond the commonly used mean-field approximation. To this aim we include neutrino-antineutrino 
pairing correlations to the two-body density matrix. The inclusion of these new contributions leads 
to an extended evolution equation for the normal neutrino density and to an equation for the abnor- 
mal one involving the pairing mean-field. We discuss the possible impact of neutrino-antineutrino 
correlations on neutrino flavour conversion in the astrophysical and cosmological environments, and 
possibly upon the supernova dynamics. Our results can be easily generalized to an arbitrary number 
of neutrino families. 
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I. INTRODUCTION 

Oscillations between quantum states is a widespread phe- 
nomenon, appearing in different physical contexts, like 
the Rabi oscillations in optics, or the Kg — Kq and neu- 
trino oscillations in particle physics. The propagation 
in a medium can produce resonant conversion between 
quantum states, such as when neutrinos change their 
flavour while travelling in a star, or when photons mod- 
ify their polarization in a bircfringcnt medium. The ba- 
sic oscillation phenomenon can be modified sometimes in 
surprising ways, if the interaction with a medium intro- 
duces complexity. The study of how neutrinos change 
their flavour in stellar environments and in the early uni- 
verse has uncovered numerous such examples. Complex- 
ity arises from the non-linear and the many-body char- 
acter of the problem. 

Neutrinos are elementary particles having non-zero 
mixings, as first conjectured by Pontecorvo |l| and dis- 
covered in 1998 by Super-Kamiokande Numerous 
experiments have contributed to the measurement of the 
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neutrino mixing angles and squared-mass differences that 
shape the way they change their flavour while travel- 
ling In particular, precisely known are the neutrino 
mixing angles of the Maki-Nakagawa-Sakata-Pontecorvo 
(MNSP) matrix that relates the interaction (flavour) to 
the propagation (mass) eigenstate basis the squared- 
mass difference values and one sign, while the other sign 
remains unknown (the hierarchy problem). Addressing 
the question of the value of the (Dirac or Majorana) CP 
violating phases is one of the major future goals, jointly 
with the determination of the neutrino absolute mass and 
of the neutrino (Dirac or Majorana) nature [f|. 

Astrophysical and cosmological environments produce 
copious amounts of neutrinos. Therefore, a precise 
knowledge of neutrino flavour conversion in media is re- 
quired, to assess e.g. the neutrino impact on the su- 
pernova dynamics and on (stellar or primordial) nucle- 
osynthesis processes, to interpret the signal associated 
with solar neutrinos, to predict the one produced by 
core-collapse supernovae, or to understand how neutri- 
nos change their flavour while traversing the Earth. It 
is experimentally established that the origin of the so- 
lar neutrino deficit is a resonant flavour conversion in- 
duced by neutrinos interacting with matter. This is the 
well known Mikheev-Smirnov-Wolfcnstcin (MSW) effect 
0, [?J . While this is the reference phenomenon to under- 
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stand how neutrinos change their flavour in media, var- 
ious other phenomena impacting flavour conversion can 
occur, depending on the specific environment under con- 
sideration. For example, Pantaleone has first pointed out 
the presence of a non-linear refractive index due to the 
neutrino interaction with other neutrinos whenever the 
neutrino number density is large Q. This is relevant for 
neutrino evolution in the early universe, in core-collapse 
supernovae, and for low energy neutrinos produced in 
accretion-disk black-hole scenarios. Indeed, simulations 
implementing the neutrino-neutrino interaction show the 
emergence of new phenomena as pointed out by Samuel 
Q , that can be interpreted as the synchronization of ef- 
fective spins (loj . a flavour (ll| or gyroscopic pendulum 

12 1, or a magnetic resonance phenomenon [13[ (see e.g. 

14| for a review). Features associated with the explo- 
sion dynamics of a core-collapse supernova, such as the 
presence of shock-waves and of turbulence, produce new 
interference phenomena, like multiple MSW resonances 
fl5| and eventually depolarization [l|| Qj} . The spectral 
and location changes of the supernova neutrino fluxes im- 
pact supernova observations and are currently being in- 
vestigated. In the cosmological context, neutrino flavour 
conversion is also important at the epoch of Big-bang nu- 
cleosynthesis (see e.g. [H, [lj| for a review). Numerous 
works have investigated for example the effects of mixings 
among active flavours [20l - [22l, w ith a possible non-zero 
leptonic CP violating phase |23l . |24| or between active 
and sterile neutrinos [24l427| on the primordial element 
abundance(s). 

Numerous works have aimed at formulating theoreti- 
cally the equations of motion that describe the evolution 
of particles with mixings in a medium. The MSW effect 
is usually accounted for, by using an effective Hamilto- 
nian that is linear in the weak coupling constant, and 
depends upon the matter number density Q. In Q 
a generalization of the neutrino evolution equations is 
made in a similar way, to implement interaction of neu- 
trinos with themselves. Pantaleone already emphasizes 
the complexity inherent with the non-linearity and the 
many-body character of the problem. Ref.Q has given 
a mean-field equation including such interaction terms. 
Refs.[28| and [29[ have derived neutrino evolution equa- 
tions beyond the mean-field approximation and including 
two-body collision terms with the "molecular chaos" as- 
sumption that neglects the building up of correlations in 
the collision term. Such Boltzmann equations are formu- 
lated in terms of matrices of neutrino densities. Ref.[3~i| 
has also derived the neutrino Boltzmann evolution equa- 
tion using first quantization, having in mind the case of 
the early universe. Ref. [32T ] has generalized the equations 
of [29| to the three flavour case (without the collision 
terms) and made the angular dependence more explicit. 



Many-body aspects, and the possible breakdown of the 
one-body description, have been discussed in Refs.(33l- 
[35j . using the spin-spin analogy in simplified models. Li- 
ouville equations for neutrino distribution matrices are 
derived in [36) . An algebraic approach to the neutrino 
propagation in media is given in [37| and the evolution of 
the many-body problem is formulated as a coherent-state 
path integral. This allows one, in particular, to calculate 
corrections to the mean-field e quat ions as a determinant 
coming from the path integral (37|. The algebraic based 
formulation is further employed in [38j . where the neu- 
trino Hamiltonian of the many-body system (with mix- 
ings and the neutrino-neutrino contribution but without 
the matter term) is put in connection with the (reduced) 
BCS pairing Hamiltonian describing superconductivity. 
It is pointed out that the corresponding constants of mo- 
tion show the exact solvability of the problem. 

The present work formulates the problem of the neu- 
trino evolution, in terms of reduced density matri- 
ces, using the Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy This replaces the Liou- 

ville Von-Neumann equation for the many-body density 
matrix describing the full many-body problem, by an 
(unclosed) set of coupled intcgro-differcntial equations 
for the reduced density matrices. Using this theoreti- 
cal framework, we derive the equations of motion for a 
system of neutrinos, travelling in a medium made of or- 
dinary matter, and interacting with each other. First we 
show that truncating the hierarchy at lowest order pro- 
duces the mean-field equations for the one-body density 
matrix, commonly used in the literature, to implement 
the neutrino interaction with matter and with neutri- 
nos. Next, we focus on the neutrino evolution descrip- 
tion beyond the mean-field approximation and include for 
the first time neutrino-antineutrino pairing correlations 
to the two-body density matrix. Such bilinear products 
have been neglected so far, based upon the argument that 
their expectation value over free states typically oscillate 
fast around zero. However, since the neutrino evolution 
equations are often non-linear, it is worthwhile to investi- 
gate their possible impact on neutrino flavour evolution 
in a medium. We show that the inclusion of neutrino- 
antineutrino pairing correlations leads to extended time- 
dependent mean-field equations both for the normal and 
for an abnormal neutrino density matrix. Finally we 
conclude by discussing the possible implications of these 
contributions for neutrino flavour conversion in the as- 
trophysical and cosmological environments. 

The manuscript is organised as follows. Section II 
presents the theoretical framework of the BBGKY hi- 
erarchy, its lowest order truncation that furnishes the 
mean-field approximation and the evolution equation for 
the two-body correlation function. In Section III we red- 



3 



erive the mean-field neutrino evolution equations includ- 
ing both the coupling to matter (the MSW contribution) 
and to neutrino themselves (the vv interaction term). 
We discuss the relationship with the equations commonly 
used in the literature. Section IV introduces the con- 
tribution from vv pairing correlations to the two-body 
correlation function. Our extended mean-field equations 
for the normal and abnormal neutrino density matrices, 
involving the normal and pairing mean-fields, are pre- 
sented. Section V includes a discussion and our conclu- 



II. THE THEORETICAL FRAMEWORK 



. dpi... s 

' dt 



tr s+1 [v£i' S) ,Pi... s +i] (4) 



wher^l 



Pi...s = (at . ..a\ai ...a s ) (5) 

denotes the s-body density matrix components^. The 
a s and a\ correspond to the particle annihilation and 
creation operators for a particle in the quantum state s, 
respectively. In particular, the 1-body and 2-body matrix 
elements components are: 



Pi 



(a{ai) 



(6) 



A. The Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy 

In numerous contexts one is interested in determining 
the dynamics of a system made up of N interacting par- 
ticles, such a gas of weakly interacting neutrinos, or an 
ensemble of strongly interacting nucleons in a nucleus or 
in a collision among nuclei. The Hamiltonian for such 
a system of N-particles, interacting through a two-body 
interaction, reads 



(1) 



k<k' 



comprising the one-body Hq kinetic and the two-body 
V interaction terms. The k, k' indices run over single- 
particle quantum states, identified by single-particle 
properties like momentum, flavour, helicity, isospin, etc... 
The system's evolution is determined by solving the 
Schrodingcr equation for the many-body quantum state 
\tp{t)} in case of a pure state or, more generally, the Liou- 
ville Von-Neumann equation for the many-body density 
matrix D: 



.dD 

'~dt 



[H,D] 



(2) 



with H given by Eq.([T]) (here we take U = c = 1). In the 
BBGKY hicrarchjo theoretical framework [33 - l42l |. one 
introduces the s-reduced density matrix p 1 " s defined as: 



AH 



P 



(N-s)\ 



rtr 



S+1...N 



D 



(3) 



tr s+ \ indicating that we are tracing over the s+1 particle, 
and replaces Eq.@ for D, by an unclosed equation for 



1 See e.g. Ref.[43 



P12 = {a\a\aia 2 ) (7) 

In Eq.(|l]) is the Hamiltonian of the sub-system of s 

interacting particles, while v}l'{ s ^ = J2k V{ k i s + with 
k=l. . . s. This equation is unclosed since the p\... s evo- 
lution is coupled to the one of the (s+l)-reduced density 
Pi...s+i, via the two-body interaction. The BBGKY hi- 
erarchy ((3]) can easily be deduced by applying successive 
traces to Eq.([2|) and using the property: 



1 



Pi. 



N 



-tfs+lP\...s+l 



{N-s)\ 



rtr 



s+1,... TV 



D (8) 



More explicitly Eq. (j4|) can be written as a hierarchy of 
equations of motion for the 1-body p\ to the s-reduced 
pi...s density matrix components^: 

' ipi = [Ho(l),pi}+tr 2 [V(l,2),p 12 } 
Wi2 = [H (l)+H (2) + V(l,2),p 12 ] 
+tr 3 [V{l,3) + V(2,3),p 123 ] 

'/<:...- - EUffoM + Ei^i^MVi....] 



Efc=i^+i[ v '( fc ' s + 1 )'Pi. 



.S+lJ 



(9) 

Solving Eq.® is completely equivalent to determining 
the exact evolution for D Eq.®. The advantage of the 
BBGKY framework is that it furnishes a hierarchy of 
evolution equations for the reduced density matrices of 
increasing order, so that one can test different approx- 
imations, by going at a higher truncation level in the 
hierarchy. 



2 Note that we do not write explicitly "hat" over the creation 
and annihilation operators through the whole manuscript, not 
to overburden the text. 

3 Note that we denote with p 1, B the operators, while we indicate 
the density matrix components with pi... B or p(l . . . s). 

4 Note that from now on we will denote d/dt with a dot. 
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1. The mean-field approximation for the evolution 
equations 

Let us consider the first equation of the BBGKY hierar- 
chy: 



if,! = [H (l), Pl }+tr 2 [V(l,2), Pl , 



(10) 



One can separate the correlated from the uncorrelatecH 
contribution of the two-body density matri:>|5 



Pl2 = P\p2 + Cl2 



(11) 



where en is the two-body correlation function. If one 
deals with identical fermionic particles the uncorrelatcd 
contribution p\p2 has to be replaced by p\Pi(\ — P12), 
with Pi 2 being the operator that exchanges particle 1 
with particle 2, to properly account for the antisym- 
metrization. Inserting Eq. (fTT|) in Eq. (fl"U)) . one obtains: 

i Pl = [P (l), Pi] + tr 2 [V{l, 2), p lP2 ] + tr 2 [V(l, 2),c 12 ] 

(12) 

Here, since no approximation is made, the dynamical 
equation for the one-body density is exact. 

Now, neglecting the correlated contribution to the two- 
body correlation function, one gets: 

ipi = [P (l),pi]+<r 2 [1/(l,2),pip 2 ] (13) 

or, equivalently: 

ipi = [hi(p),Pi] (14) 

with hi{p) = ff (l) +Ti(p) and r x (p) = tr 2 (V(l, 2)p 2 ) 
being the mean- field acting on particle 1. This is the 
so-called mean-field approximation. 

Writing such an equation more explicitly, it reads: 



if>i,ij ~ [Hq(1) + Ti(p), pi]j 







with 



Fi.tj(p) 



^ ' v (im,jn)P2,nm 



(15) 



(16) 



The mean-field potential is built up from a complete 
set of one-body density matrix components for parti- 
cle 2 p2,mn, each contributing with the matrix elemenlQ 



5 We will also use "linked" and "unlinked" to denote the correlated 
and uncorrelated contributions, respectively. 

6 Note that here p2 = Pi (2) is the one-body density matrix asso- 
ciated with particle 2. 

7 Note that in case of identical particles the matrix elements are 
antisymmetriscd, i.e. V(im,jn) = (im\V\2\jn) — (im\ V12 \nj). 



V(im,j n ) = (im\V 12 \jn), with jn (im) incoming (outgo- 
ing) single-particle states. From Eq. (fT6|) one can see the 
dependence of our mean-field on the one-body density 
associated with particle 2, while in some cases the inter- 
action itself might also have an explicit dependence on 
p. To solve Eq. fTS"]) one has to assign the state of the 
many-body system at initial time, which can be either 
a correlated state, or a product of independent single- 
particle states. In the latter case, the condition inherent 
to Eq.tfTB")). i.e. C12 = 0, ensures that it stays as such at 
any time. 

It is worthwhile to mention that a first-order evolution 
equation for the one-body density matrix associated with 
a given D can also be obtained by applying the Ehrcnfest 
theorem: 



= (i a ] a iiH}) 



(17) 



with pi t ij = (a^cii) and H given by Eq.([T]) (in second 
quantization). In particular, Eq. (fT5]) is recovered when 
neglecting the correlated contribution to the two-body 
density matrix [44[ . 



nd the mean-field approximation 



2. 



Our main goal will be to discuss contributions beyond the 
mean-field approximation given by Eq. (|14[) , to the evolu- 
tion equations for a system of relativistic neutrinos that 
interact among themselves and with matter. To this aim 
a useful reformulation of Eq.Q is given by a hierarchy of 
evolution equations for the correlation functions (details 
of the demonstration can be found in Ref . [45| ) . where 
only linked terms are shown to remain. Such a reformu- 
lation has the advantage that higher-order contributions 
are expected to decrease with increasing rank 45(. In 
this context, the mean-field equation (|14p is unchanged; 
while one gets for the two-body correlation function [45| : 

zci 2 = [hi(p) + h 2 (p),cn] 

+ (1 - pi)(l - P2 )V(1, 2) PlP2 (l - P i2 ) 

-(1 - Pi 2 )pip 2 F(l, 2)(1 - pi)(l - p 2 ) 

+ (1 - Pi - P2)V{1, 2)cia - c 12 V{\, 2)(1 - pi - P2 ) 

+tr 3 [V(l, 3), (1 - Pi 3 )pic 23 (l - P12)] 

+tr 3 [V{2, 3), (1 - P 23 )p2Ci 3 (l - P12)] 

(18) 

with h\(p) (h 2 (p)) the mean-fields acting on particles 1 
(2) respectively and Pi 3 (P 23 ) is the operator that ex- 
change particle 1 (2) with 3 (see Appendix A for an ex- 
plicit formulation of Eg- fTS]) ). Such an equation contains 
three main contributions coming from two-body interac- 
tions. Retaining the second and third term on the r.h.s. 
of Eq. lfT8|) . and making the "molecular chaos" assump- 
tion, that the build up of correlations due to collisions 
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is negligible, one gets a collision term with incoming and 
outg oing particles described by free particle states (see 
e.g. |46[). Such a term gives rise to a Boltzmann equa- 
tioioThe (1 — pi)(l — P2) factoi@ ensures the appropri- 
ate statistics (no contribution if the final single-particle 
states are already occupied). The fourth and fifth terms 
on the r.h.s. of Eq. (|18[) implements contributions from 
the correlated part of the two-body correlation function. 
It is on this term that we will focus later on. In partic- 
ular we will consider cases involving both particles and 
anti-particles. In such systems, one might have possible 
contributions from the expectation values of the product 
operators of the type a k bi and a k b\. Such bilinear prod- 
ucts include particle-antiparticle correlations that can be 
seen as correlations of the pairing type. These terms have 
been neglected so far, based on the argument their ex- 
pectation value (over free states) typically oscillate fast 
around zero, (a* (p,t)tf(j)',t)) (see e.g. [H, HI)- Fi- 
nally, the three-body terms give the contribution from 
the two-body interaction among three particles, obtained 
by tracing over the third particle. We will neglect these 
higher order correlations here. Since, in this work, wc 
focus on the inclusion of vv contributions to ci 2 , we will 
not consider the collision and the three-body terms. Our 
evolution equation for the two-body correlation function 
is: 



ic 12 = [hi(p) + h 2 (p),c 12 ] 

+ Pl - p 2 )V(l,2)c 12 



(19) 

ci 2 V(l,2)(l-pi-p2) 



B. The application to neutrinos 

We are here mainly interested in discussing the neutrino 
evolution and flavour conversion when neutrinos prop- 
agate in an astrophysical environment, or in the early 
universe. While the BBGKY hierarchy is generally em- 
ployed for systems of interacting particles without mix- 
ings, we here consider that the density matrices in Eq.Q 
are associated with mixed particles. Note that a "matrix 
of densities" , generalizing the usual occupation numbers, 



8 Note that a Boltzmann equation for a system of neutrinos and 
anti-neutrinos is derived in [28l |29| and in [3l{ . Such evolution 
equations include a collision term as required in the context of 
the early Universe. In core-collapse supernova simulations, full 
transport equations for neutrinos (but without the inclusion of 
mixings) are usually employed in the dense region where neutri- 
nos are trapped (see e.g. |48b52|1 '). 

9 This factor as well as the (1 — pi — P2) one come from the linked 
contribution of the trace term over the third particle. 



is commonly used in the literature (see e.g. [28], |29(). 
More explicitly, in the three flavours case the neutrino 
density matrix reads 



Pv 



,t 



( a u[},j a v~ r k, 



(20) 



The off-diagonal (or coherent) terms are non-zero to en- 
code the presence of neutrino mixings. The neutrino oc- 
cupation number for a given v a .i flavour state is given by 

the diagonal element of the density matrix {cv v i a^ a .i), 
with N„ a = J2i( a t i a v a ,i) the total occupation number. 
The definition in Eq. ([2"0")) can easily be extended to the 
case of an arbitrary neutrino families, to account for the 
presence of both sterile and active neutrinos. 

Since the systems we are interested in involve both par- 
ticles and anti-particles, one introduces a density matrix 
p v , in a way analogous to Eq. ([2"0"]) . but replacing the par- 
ticle operators a) , a with the anti-particle , b ones 2^] . 
For the sake of clarity concerning the convention used in 
the present work, we give its explicit expression: 



f(bl a 
\( b l 



> (&LA-,<> 
(4 



(bl^jb^.k) 



<<, 
(bl 



,k;bvp,j) 



(21) 



To implement the anti-particle degrees of freedom, we 
consider the usual expansion of the neutrino fields in the 
Schrodinger picture 



E 

h 



d 3 p 



(2n) 3 2E p 



[a(p, fyuxhe^ 



(22) 



The summation is over the helicity h states, the inte- 
gration over the momenta p of the (anti-)particles and 
u p,h (vp,h) are the usual Dirac spinors, Ep being neu- 
trino energy. Note that in such an expression the particle 
creation and annihilation operators are associated with 
states of a given mass for which the usual anticommuta- 
tion rules are properly defined. In the present manuscript 
we emplo£3 

{a(p,h),a?(p',ti)} = (2Tr) 3 2E p 6 3 (p-p')5 hh , (23) 



and 



{a(f,h),a(p',h')} = 0. 



(24) 



10 Similar relations hold for the anti-particle annihilation b(p, h) 
and creation b'(j>', h') operators. 
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The single-particle states associated with neutrino mass 
eigenstates are 



,1) 



(25) 



with |) being the vacuum state defined by a m \) = 0. 
The flavour eigenstates are related to the mass eigen- 
states through \v a ) = J2i (i and a correspond to 
an arbitrary number of neutrino families) , where U is the 
Maki-Nakagawa-Sakata-Pontccorvo (MNSP) unitary ma- 
trix. In three flavours, the MNSP matrix depends upon 
three neutrino mixing angles that have been measured, 
one Dirac and two Majorana unknown phases @ . Since 
the appropriate anti-commutation rules only hold for the 
mass eigenstates operators, Eas. (|23W2"3|) require special 
attention (see e.g. Ref.[36| for a discussion). This sub- 
tlety is sometimes avoided by writing the fields Eq.(|22p 
for massless neutrinos, as done e.g.[2{§. As discussed for 
example in Refs.[36|, [53|, it is not possible to rigorously 
build up a Fock space for flavour states since the neu- 
trino flavour creation and annihilation operators do not 
satisfy the canonical anticommutation rules Eas. (|23M2"3]) , 
Indeed, as shown in Ref . [531] . one can define an approxi- 
mate Fock space by introducing neutrino "weak-particle 
states" that depend upon the specific weak process un- 
der consideration. However, in the limit of rclativistic 
neutrinos, the anticommutation relations (|23H24[) also ap- 
proximately hold for flavour states. Since this is a good 
approximation for our cases of interest (solar and super- 
nova neutrinos, cosmological neutrinos at the epoch of 
big-bang nucleosynthesis), we make the assumption that 
an approximate Fock space for our flavour states can be 
built. It is worthwhile to mention that extending our 
results without making this approximation does not in- 
troduce extra conceptual difficulties. In much the same 
way as done in Ref. 53] for the muon decay case, but in 
the mean-fields expressions, one should retain an explicit 
dependence on the mixing matrix elements in the calcu- 
lation of the interaction matrix elements, depending on 
the specific weak process under consideration^ . 

Finally, it is worthwhile to remind that, while the 
BBGKY theoretical framework is based upon the density 
matrix formalism, the evolution equations §§§ can be for- 
mulated as a hierarchy for many-body Greens' functions 
as done in Ref. [47j . In particular, in the equal time limit, 
when only linked contributions are retained, these two 
formalisms are completely equivalent. 



III. NEUTRINOS EVOLVING IN A MEDIUM IN 
THE MEAN-FIELD APPROXIMATION 

We now take the example of neutrinos interacting with 
the electrons, protons and neutrons composing a medium 
to show how the formalism presented in Section II can be 
used to rigorously derive well known neutrino evolution 
equations. We just sketch the main lines of the deriva- 
tiorF^I. The assumption that is usually made is that our 
system of neutrinos interacting with a medium can be 
described at lowest order as a system of independent par- 
ticles, so that implicitly the problem reduces to following 
the evolution of a single-particle at a time and calculating 
the evolution of the associated one-body density matrix. 

Let us consider the BBGKY hierarchy truncated in the 
mean- field approximation given by Eq. (|14[) . The Hamil- 
tonian Eq.([T]) for our case of interest in the flavour basis 
reads: 



H f = UH m U ] + H„ 



(26) 



where U is the MNSP unitary matrix relating the neu- 
trino flavour basis to the mass eigenstate basis (with 
eigenenergies Ei). The first H m = diag(Ei) contribution 
is the propagation term, while the second one Hi nt cor- 
responds to the two-body interaction between a neutrino 
and another particle. 



A. The mean-field associated with neutrino 
interaction with matter (MSW case) 

To follow the one-body density matrix evolution given 
by Eq. p"?]) . one needs to determine the mean- field (fll))) 
created by the background particles and acting on the 
"test" neutrino. The interaction term ff„j corresponds 
to the charged- or neutral-current Hamiltonian describ- 
ing the neutrino interaction with the medium. We take 
the example of the charged-currcnt interaction on elec- 
trons, where Hi nt is given bj0 



cc 



Gf 

7f 



d 3 x{cj) eltl {l - 7 5 )^J[0^7 M (1 - 7s) 



(27) 

where Gp is the Fermi coupling constant and with 
the fields <f> given by Eq. (|22|) . This requires calculat- 
ing the matrix elements for neutrino-electron scattering 



11 To implement this correction, one should keep in mind that we 
deal with the process amplitudes, and not amplitude squares as 




A different derivation of the mean-field equations accounting for 
the neutrino interaction with a medium is given e.g. in [f|. 
The effective low energy approximation is sufficient for the ap- 
plications envisaged. 
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FIG. 1. Neutrino interaction with electrons and the corre- 
sponding mean-field associated with the electron background. 
In the evolution equations, the mean-field acting on a single 
neutrino state is build up from the summation of the electron 
single-particle states (see Eqs.JTBJ, (|28[l and (1321) ). 



v & ie m \Hcc\ l 'ej&n) with the interaction Hamil- 



tonian (f2"T)) . One gets 



d 3 p 



2tt) 3 2E p J (2n) 3 2E p , 



(2ir) 3 5 3 (p + k -p' -k') 

[U e (p, h e )^(l - 75)«e(p', ^e)] 

(ot^hJoeCT.ft')) (28) 

The summation over the m, n single-particle states in 
Eq. (jl6p becomes here a sum over the electron helicity 
states and an integration over momenta (Figure [1} . 

The key quantity to define is the expectation value 
of the aUp, h)a e (p' , h') operator, over the electron back- 
grounco characteristic of the considered mediunF^l : 



Pp'h',ph = (<4(£ h)a e (p', h')) 



(29) 



The assumption of a homogeneous and unpolarized 
medium corresponds to: 

Pp>h-,ph = (2ir) 3 2E p S hh ,S 3 (p~ p')p} (30) 

Eq. (|29l) constitutes for example a good approximation for 
the case of the Sun. Using Eqs. (|29H3"0|) the total electron 
number is 



N K = 2V 



d 3 p 

W) 



3 pp 



(31) 



14 Here we make the assumption that the contributions from the 
neutrino and the electron backgrounds can be factorized in D, 
i.e. that modifications of the electron background coming from 
the interaction with neutrinos can be neglected. 

15 Note that in this and in the following section, we denote the 
single-particle density pi kl as Phi to simplify notations. If, a 
particle is identified by helicity and momentum, then p^i reads 
Pfih p'h' f° r an °l cc tron background, for example. 



where the factor of 2 comes from the summation over 
the electron helicity states, and V is the quantization 
volume. In Eq. (|28p . by tracing over the spinors, imple- 
menting that neutrinos have only one possible helicity 
state, and assuming the electron background is homoge- 
neous, unpolarized and isotropic, one gets for the mean- 
field Eq.CS 



(32) 



T Ua (p e ) = (2ir) 3 2E k 6 3 (k - k')V2G F n e 



where n e = N e /V is the electron number density. In 
Eq. ([3"2")l the 5-function ensures that the momentum of 
the neutrino propagating in the electron medium is un- 
changed, as a consequence of homogeneity. The prefactor 
{2ir) 3 2Ek are present because of the chosen normalization 
of the (anti-)particle anticommutation relations Eqs.(|23I 
[24]) . Our result (|32|) is the mean-field contribution to 
the neutrino Hamiltonian, corresponding to neutrino in- 
teraction with matter, that is linear in the weak cou- 
pling constant and depends upon the number densities 
of the particles composing the medium. We find the well 
known low energy effective Hamiltonian (see e.g. Q) 
that can give rise, if neutrinos propagate adiabatically 
in a medium, to a resonanlF^I flavour conversion - the 
Mikheev-Smirnov-Wolfenstein effect @, 0] ■ 

The procedure just outlined can be applied to the 
neutral-current v e scattering on electrons, to the charged- 
and neutral-current v e scattering on electrons, of v e or 
v e on positrons, as well as neutral-current scattering on 
protons and neutrons, giving the expected results. In 
particular, it is immediate to show that, when dealing 
with anti-particles, the associated mean-field T depends 
upon matrix elements that involve the w and b operators 
(instead of the particle operators), introducing a minus 
sign. For example, for the case of v e evolving in an elec- 
tron background, one obtains the expected result 

^uAPe) = -(2ir) 3 2E k 5 3 (k - k')V2G F n e . (33) 



B. The mean- field associated with neutrino 
self-interactions 

The second case we are going to consider, in the mean- 
field approximation, is when the (anti-)ncutrino is evolv- 
ing in a background of v and v. This case is of interest 
since recent studies have shown the important role of 



16 The occurence of such a resonant phenomenon depends upon 
neutrino properties (energies, mixing angles, squared-mass dif- 
ferences value and sign) and the specific matter number density 
profile for the system under consideration. 
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the neutrino-neutrino interaction for instance in a core- 
collapse supernova, and the variety of new flavour con- 
version phenomena that can arise, when implementing 
this contribution in the neutrino Hamiltonian (see e.g. 
fl4| for a review). A scries of works have discussed the 
neutrino evolution equation in presence of such terms 
B S HHH [H IH|. Note that the role of such terms 
was first pointed out in the context of the early universe 
Here we just sketch the derivation of these equa- 
tions, following the same procedure as for the MSW case 
(more details are given in Appendix B). 

The neutrino-neutrino mean-field depends this time on 
the effective low energy neutral-current Hamiltonian: 



Hnc = 



d 3 x[4>^^(l - 75)^J[0^7 M (1 - 75) 



with v y = v e ,Vn or v r 



mean-field T 

corresponding anti-ncutrinos) : 

G F f d 3 p 



(p.) 



(34) 

From Eq. ([T5)) one gets for the 
V v 

d 3 p' 



v a ,vc<{Pv) with v a,vp = v e , v^-, v T (or the 



2^2) (2tt) 3 2E p J {2tt) 3 2E p , 
(2n) 3 S 3 {p + k-p' -k') 
[u Vff (k, h^)^(l - 75)^ (k\ h' a )} 
[u„ a {p,h a )^{l - 7 5 )u^ if', tip)] 
{at a (p,h a )a Ufl (p',tip)) 
and requiring the expectation value : 



(35) 



Pp'h"ph = ( a L (P> h a )a vp (p\ tip)) (36) 
correspond to a homogeneous and unpolarized system: 



Pp'h',ph 



{2Tr) 3 2E p 6 hh ,6 3 (p 



P )Pp 



(37) 



The key difference with the case of the electron back- 
ground is that the quantity PpifJ x p l has diagonal and 
off-diagonal terms (Figure [2]). As a consequence the 
mean-field acting on a v e for example has both diago- 
nal T VajVa (p v ) and off-diagonal r * Va , V/} (p v ) contributions. 
One recognizes in the mean-field contributions that the 
diagonal one is associated with forward scattering with- 
out flavour exchange, while the off-diagonal contribution 
is associated with forward scattering "with flavour ex- 
change" (the Pantaleone diagram Q). We emphasize 
that in our derivation we do not directly start with an 
effective Hamiltonian having these terms, as done e.g. in 
Ref.[13|. Here these contributions naturally arise when 
considering the weak interaction term, using the density 
matrix Eq. (|2T)]) and calculating T (or, in other words, 
closing the loop). By implementing that neutrinos have 
one helicity state only, tracing over the spinors, in case 




v0) 




u a (k) 



v a {kf) 



u p (k) 



FIG. 2. The figure shows, in a pictorial way, the interaction 
terms and the corresponding mean-fields Eq. (|38|) associated 
with the neutrino-neutrino interaction. The mean-field act- 
ing on a single neutrino state is build up from the summation 
of the single-particle states that make up the background. 
The two contributions correspond to the diagonal part of the 
mean-field V VaiVa , arising from the usual scattering terms (up- 
per figures), and the off-diagonal part F Va:V/3 Eq. (|36[) (lower 
figures) associated with Pantaleone off-diagonal refractive in- 
dex. The diagonal contribution to the mean-field T„ atUa has 
an extra term coming from neutrinos of the same flavor a 
running in the loop. Both the diagonal and the off-diagonal 
mean-field receives similar contributions from anti-neutrinos, 
instead of neutrinos, in the loop Eq. (|39p . 



the neutrino background is homogeneous and anisotropic, 
one gets for the off-diagonal contribution the following 
expression: 



= (2ir) 3 2E k 5 3 (k - k') 
d 3 p 

W) 



V2GI 



with p = p/\p\ and k = k/\k\. A similar expression holds 
for the diagonal one if Vp = v ai except for an extra factor 
of 2. In the case of an isotropic medium the angular term 

cosQ^~ = k -p averages out, and one recovers the same re- 
sult as Eq. (pl2"]) . If the background includes anti- neutrinos 
as well, one needs to add the contribution from p* , but 
with a minus sign. By adding up the two contributions 
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the total neutrino-neutrino mean-field readeFI: 

T{p v ,p v ) = {2Tr) 3 2E k 5 3 {k- k')V2G F (39) 

where here p Vat p (p* ^) stand for the density matrix 

Eq. (f2T)| (Eq. ([2"TjO . v a refers to a neutrino that is ini- 
tially born in the a flavour. In fact, one has to sum over 
all neutrino flavours present in the system. Note that in 
Eq. (f39| the trace term, tr{p Va $ — p*^ -), has been sub- 
tracted. 



homogeneity. In case the quantum state at initial time is 
an independent particle state, the one-body density ma- 
trices (|20M21[) can be replaced by single-particle one-body 
densities. In the neutrino case, the diagonal elements pi 
directly give the neutrino survival probabilities \vi\ 2 Vi 
being the neutrino amplitude for flavour i, while the off- 
diagonal ones pij are the mixing terms v^i/j (and simi- 
larly for anti-neutrinos). We conclude here by empha- 
sizing that the evolution equatio ns (14011431) we find are in 
agreement with those of Refs.j29ll30l.l37j]. commonly used 
in the investigation of solar, of supernova neutrinos and 
of the low energy neutrinos in accretion-disks black-hole 
scenarios. 



C. Neutrino evolution equations in the mean-field 
approximation 

With the results of Eqs.([22|) and (EH, the mean-field 
equations (|15M16[) for the density matrix (|20p describing 
the neutrino evolution in a medium becomes explicitly: 

ip=[h(p),p] (40) 

with 

h(p) = UH m U^ + H mat (41) 
+V2G F £ J ^(P^, P ~ P^, P ) (l - P- k) 

where H mat = diag(V2GFn e , 0,0). A similar equation 
holds if an anti-neutrino is travelling instead of a neu- 
trino: 

ip=[h(p),p] (42) 

witl0 

h{p) = U*H m U T - H mat (43) 

-V2G F J2 J (^^Lp " i 1 ~ P' *) 

Note that in Eq . (|40H4^|) p and p have two indices in 
flavour, but only one in momentum as a consequence of 



17 Note that the contribution coming for the anti-neutrino density 
matrix is p* if one employs the definition pij = (b^bi) (see Ap- 
pendix B). If one defines such a quantity as pij = {b\bj) then the 
contribution appearing in the evolution equations depends upon 
p, instead of p* . 

18 Note that, according to our definition for p Eq J21ll , anti- 
neutrinos do not transform the same way as neutrinos under 
U. 



IV. NEUTRINO- ANTINEUTRINO PAIRING 
CORRELATIONS 

A. The extended neutrino dynamical equations 

Our main goal is to obtain the evolution equations, be- 
yond the mean-field approximation, for a system of neu- 
trinos and anti-neutrinos evolving in an environment, 
taking into account possible vv pairing correlations. The 
latter correspond to the following linked contribution to 
the two-body density matrix Eq. (|ll[) : 

c a /3, a 'p' « (a a 'bp,){bpa a ) = K* a , p ,K al3 (44) 

where the quantities K a p and k*/^/ are called the ab- 
normal densities. Imposing that the pair products con- 
serve individual lepton numbers, only product of oper- 
ators associated with the same flavour, such as a Va b Va , 
and the corresponding hcrmitian conjugates, are admit- 
ted. Since neutrinos have mixings, one can have con- 
tributions from products involving neutrino-antineutrino 
pairs with different flavours. More generally, e.g. if to- 
tal lepton-number is not conserved one could have cor- 
relations associated with the pair product operators like 
Q>v a Q<v a (and similarly for antineutrinos) . We note that 
the expectation value J*^l (a^a), (b^b), (ba) and {a%^) nat- 
urally appear as components of the field correlation func- 
tion (4>4>) ■ So fax, neutrino-antineutrino correlations have 
been neglecteco (see e.g. Refs.(28l. [29l. |36|). For exam- 
ple, in the formulation of Ref. [36| , they correspond to the 



Note that the spinor products are not usually modified in our 
cases of interest. 

Note that, in the context of baryogencsis via leptogenesis, 
the authors of Ref. |55j have emphasized the role of neutrino- 
antineutrino correlations, in a quantum field theory approach, 
including the Boltzmann collision term but within a simplified 
neutrino model. 
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rapidly oscillating cross terms between the positive- and 
negative- frequency parts of the quantum density function 
iTlj 1 — {Nv l i {y)V 7 p{z)), where the v{y) and i>(z) are the 
neutrino and anti-neutrino quantum field operators. 

Two perspectives are possible to investigate the impact 
of the neutrino-antineutrino correlations on the neutrino 
evolution. The first is to assume that such correlations 
are non-zero at initial time, e.g. at the neutrino sphere in 
a supernova, from previous interactions among neutrinos 
in the dense supernova region, where they are trapped. 
The second possibility is to study whether such terms 
can be dynamically produced through the interactions. 
However this is more demanding, since one has to retain 
the collision term in the two-body correlation function 
Eq. (fT8"| as well. In this manuscript we adopt the first 
perspective. We will see below that, according to our ex- 
tended evolution equations if their contribution is zero at 
initial time it is zero at all times. Therefore, in this case, 
the mean- field approximation given by Eqs. (|40ti43[) is cor- 
rect (if for the considered system, the collision term can 
also be neglected in Eq.(|lgp). This is the approximation 
that is usually implicitly made. 

Let us now discuss how the neutrino evolution equation 
Eqs. (|40j|45)) have to be extended to implement pairing 
correlations between v and v. In this case, the evolution 
of the system is determined by using the first two equa- 
tions of the BBGKY hierarchy Eq.©. By substituting 
Eq.([4"4"|) in Eq. (TT9")) one obtains the evolution equation 
for the abnormal density (and its complex conjugate), 
while from the first equation of the BBGKY hierarchy 
and Eq. pTj) . we obtain the evolution equation for the 
normal density. We finally get the extended evolution 
equations: 

' ipij(l) = Kl),p(l)]«+Em( A imKj m -«imA| m ) 

( iflw(2) - [h(2),p(2)} kl +j: m (A mk K* ml -K mk A* ml ) 
iKik = Emt^mf 1 )^' + h km {2)n im ) + A ik 

{ -E m (Pim(l)A mk \ Pkrn 

(2)A iTO ) 

(45) 

where here the indices i,j stand for v a ,vp\ k,l for v a , i/g 
with a, (3 that vary over the different electron, muon and 
tau flavour states. For the sake of clarity, in Egs.(|45p we 
show explicitly the dependence on particle 1 and particle 
2 of the quantities. Obviously, in our extended equations 
one has to determine the evolution of the normal density 
associated with a neutrino (particle 1) and an antineu- 
trino (particle 2) to determine the neutrino-antineutrino 
pair evolution in k and n*. We emphasize that, the new 
evolution equations Eq. (|45|) can also be derived by using 
the Ehrenfest theorem and determining the first-order 
evolution equations for the abnormal density: 

ik ik = ([b k ai,H]) (46) 

and for the normal one using Eq. fTT]) . 




FIG. 3. The figure shows, in a pictorial way, the interaction 
of a neutrino with an anti-neutrino and the corresponding ab- 
normal mean-field A* Eq. (|48|) which is build up from the sum 
of single-particle neutrino-antineutrino states. Its expression 
is analogous to the one of the normal mean-field by replacing 
the normal density by the abnormal one re and by summing 
over the final (initial in A) states. The contributions shown 
correspond to the diagonal component A* q v Eq. (|56l) (upper 
figures), and to the off-diagonal one A% atB Eq. (|52[) (lower 
figures). Note that the diagonal abnormal mean-field also re- 
ceives a contribution from a term where neutrinos of the same 
flavor a run in the loop. 

The expression for the abnormal mean-field is: 

A* = ^2v( ikjjl )Kji (47) 
jl 

and of its complex conjugate 

A *ik = J2 V Ul,ik)Kjl (48) 
jl 

One can see that such expressions are analogous to the 
one for the mean-field T Eq. tfTB"]) , but with the abnor- 
mal density k replacing the normal density p, and by 
summing over the initial (or final) single-particle states 
instead of over a final and initial single-particle state. 

We note that Eqs. (|45|) can be cast in an elegant and 
compact matrix form: 

ili=[H,1l] (49) 
where we have introduced the generalized density 1Z: 

*-Ui-v) <5o) 
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that depends upon both the normal densities for the neu- 
trinos and the anti-neutrinos and the abnormal density. 
The generalized Hamiltonian H governing the evolution 
is given by: 



H 



h A 

At -h* 



(51) 



It comprises the mean-field Hamiltonians for the neu- 
trinos and the anti-neutrinos, as well as the abnormal 
mean-field. 



B. The abnormal mean- field 

Let us now compute the expression of the abnormal 
mean-field A. First, the expression for the mean-field 
h (h*) acting on a neutrino (anti-neutrino) is given by 
Eqs. (pTTj) and To derive an explicit expression for A 
Eq.(|j7]l and A* Eq.([i5|). one needs to calculate the ma- 
trix element v^j t ki) associated with the neutral-current 
interaction of a neutrino with an anti-neutrino. One ob- 
tains the following expression for the off-diagonal com- 
ponents of the abnormal field (see Figure [3]): 



A^«) = - 



Gi 



d 3 p 



d 3 k 



2V2J {2tt) 3 2E p J (2ir) 3 2E k 
(2ir) 3 6 3 {p + k -p' -%') 

[u Ua (k, /i q )7m( 1 ~ J5)u Ua {k', h' a )] 
[vv if, h'ph^il - 7sH (3 {p, hp)} 
(at a (k,h a )bi e (p,h )) (52) 



The explicit expression of the abnormal field Eq. ([52 
depends on the properties of the background: 



kh a ,ph@ 



(at a (k,h a )bt (p,hp)) 



(53) 



Different assumptions can be made on the expectation 
value of the correlator al a (k, h a )b\, e {p, hp) and its com- 
plex conjugate. In order to remain as general as possible, 
various options are considered in the following. We give 
the results for A* q p while similar expressions hold for 



as well as for the diagonal contributions of the 



abnormal mean-field A„ qjPq and its complex conjugate 
A* q Pq (only the final results are presented here, while 
more details on their derivation can be found in Appendix 
C). In particular, for the calculation of the latter, two 
contributions need to be added: 



{<) = 2A^ % a 3 «) + A^g, «) (54) 
2v 



E 

k,l 



* 1 a,P * -I 

{H,ij) K V<x,k,V a ,l + V (kl,lj) K Vf3,k,Vl3,ll 



(55) 



as can be seen from Eq. p%|) . The first contribution 

^■Vai,v a ■ comes from the two amplitudes for the pro- 
cess of neutrino-antineutrino scattering and annihilation 
for equal flavours (see Figure [3]) that summed give a fac- 
tor of 2 times an expression similar to the one shown in 
Eq.®: 



► (eg) 



G, 



d 3 p 



d 3 k 



2V2J (2tt) 3 2E p J {2n) 3 2E k 
(27t) 3 S 3 (p + k-p' -k') 
[u Va {k, /j q )7 m (1 - J5)uv a (k', h' a )} 
[v Va (p',ti a )^(l --y 5 )v Ua (p, h a )] 
(at a (p,h a )bt a (p,h a )) (56) 

The second contribution to A^"™^ comes from different 
flavours than a running in the loop (Fig(3J). Let us now 
perform an explicit calculation of the pairing mean-field, 
depending on the background properties. 

The requirement that the background through which 
neutrinos are traveling is homogeneous, is fulfilled, if 
we impose that the neutrino and the anti-neutrino have 
opposite momentum in the expectation values of the 
neutrino-antineutrino pair operators. This corresponds 
to: 



k ■ I ' 

and we take h a ■ 



(2ir) 3 2E k 6 3 {p 



(57) 



1 and hp = +1. Note that the 
a *) are related to the neutrino- 
antineutrino pair number density: 



quantities K~ vp 



(or 



dp v a vp* 
J2^ K p 



(58) 



in a way analogous e.g. to the electron number densities 
(see Eq. (pITj) ). From Eq. ([5^1) . by employing the procedure 
discussed in Section III. A and III.B for the normal fields, 
one obtains that the abnormal mean-field in spherical 
coordinates reads: 



^(nt) = -N(k',p')V2G 
dOdpd<p ' 



(59) 



(1 + cos 9 cos 9') cos( 



(2tt) 3 

+ i(cos 6 + cos 9') sm(<fi' — 



sin ft sin t 



sin 9p Kp 



,1/(3* 



with N(k',p') = (2n) 3 2E k ,S 3 (p' + k r ). Note that in case 
there is cylindrical symmetry only the last term remains: 



,Vf> 



N(k',p')V2G F 
d cos 9dp 
~J2^ 



sin 9 sin 9 p /c J 



(60) 
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The expression Eq. (|5i?]) gives in cartesian coordinates: Eq. ([4"5")) gives to the following set of equations. 



«) = -N{k',f)y/2G! 



d 3 p [d R (p,p') + 



(2tt) 3 



\p\\p'\\ P. 



pUp'x 



; K f> 



with 



dn(p,p') =(|pl|p'| + PzP' z )(PxP' x + PyPy) 



+ (pt+Pi)(Px+Py) 



'2\ 



and 



d l(p,p') = (PzW'\ + PZ'\P\)(P' V PX -PxP'y) 



(61) 
(62) 

(63) 



So far we have neglected the neutrino masses so that 
v and v have definite helicities. Let us now consider the 
case that e.g. h a = hp = +1 or — 1: 



k- _ 

k,p 



(2Tr) 3 2E k S ha j lfi d 3 (p + k)4 c 



(64) 



While such a contribution is possible, we expect the cor- 
responding abnormal field to be suppressed by a factor 
(m/E) 2 , E being the neutrino energy. For this specific 
derivation, since one is sensitive to the neutrino masses, 
the calculation has to be performed by working in the 
neutrino mass eigenstate basis. Following this procedure, 
one finds: 



A* ,,«) = (2n) 3 2J E%,EU 3 (p' + k') 



k' 

d 3 k mjirriA 



2y/2 



(65) 



(2n) 3 El El 



?{<!>' -<t>) (\ 



COS U%p)K j. 1 



where here the notation Ei, indicates the energy of a 
neutrino with mass ?rv and momentum k' . 



C. The three-flavor neutrino case 

We conclude by giving the expressions of the extended 
dynamical equations for three neutrino flavors. In this 
case the abnormal density is defined as: 

(K e a Ve ) (K„a Ve ) (b Vr a Ve ) 
(b Ve a v J (b Vli a Vll ) {b„ T a Vll ) | (66) 
(bu e a Vr ) (b Vll a VT ) {K T a Vr ) 

and its complex conjugate as: 

«6tJ (aiAj «t>l) 

4X> « b V « b D I (67) 
(at bt ) (at 6t ) ( a t 6 t ) 



-A 



A 



l/ a , V4 P 



(68) 

where the indices a,/3,i = e,p,,T, h (h) is given by 
Eqs. (|4T)| (Ea.(|42p) for neutrinos (antineutrinos), while 
the k v and k* are given by Eq. (|66H67)) . 

In order to determine the neutrino flavour evolution, 
one needs to assign initial conditions for the normal and 
abnormal densities, i.e. po = p(t = 0), po = p(t = 0) 
and kq = n(t = 0) for the extended mean- field equa- 
tions Eq. pgj) with (|50H5ip . In the case of core-collapse 
supernova neutrinos, such an assignment could be done 
by extracting the relevant information from realistic sim- 
ulations of the dense region where neutrinos are trapped 
that would include the relevant correlations. Otherwise, 
a guess that one can make is to assume that the system is 
approximately described at initial time by a "stationary 
state" of the extended mean- field Hamiltonian Eq. (f5~T]) ; 
while the system is driven out of such a solution at later 
times [60| . 

Finally, we would like to point out that, in the applica- 
tion of our equations to the case of supernova neutrinos, 
the conditions Eqs. ([57|) and cannot be met, strictly 
speaking. In fact, in the region outside the neutrino- 
sphere, where the density becomes low enough, all neu- 
trinos and anti-neutrinos start free streaming. However, 
in the region just before the neutrino-sphere and at its 
boundary the presence of collisions might produce non- 
zero neutrino-antineutrino correlations as we consider 
here. It is this intermediate region, that we are inter- 
ested in, where one goes from the regime where neutri- 
nos are trapped requiring a Boltzmann treatment, and 
the one where only forward scattering becomes relevant 
and the mean- field approximation given by Eqs. (|40H4"3")) 
constitutes a good approximation. The presence of v — v 
correlations in this transition region might indeed play 



21 Note that here the quantities depend upon only one momentum 
index. 
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a role from the point of view of the supernova dynam- 
ics, or eventually modify neutrino flavor conversion when 
neutrinos start free streaming. 



V. DISCUSSION AND CONCLUSIONS 

In the present work, we have derived the neutrino evo- 
lution equations in a medium composed of matter and 
eventually of neutrinos, having in mind environments 
that produce copious neutrino amounts such as the Sun, 
core-collapse supernovae, accretion disks around black 
holes or the early universe. In particular, observations 
of such neutrinos require a precise understanding of the 
corresponding neutrino number fluxes and spectra, their 
modifications through such media, and their specific sig- 
natures both in solar and supernova neutrino detectors, 
as well as on stellar, or primordial, nucleosynthesis abun- 
dances. 

In the astrophysical environments the neutrino flavour 
evolution is essentially treated in the mean-field approx- 
imation. While this is a good description in the case 
of solar neutrinos, the investigation of neutrino flavour 
evolution in core-collapse supernovae might require go- 
ing beyond. For example, in the transition between 
the dense region where neutrinos are trapped, and the 
neutrinosphere where they start free streaming, many- 
body correlations might have sizeable effects and impact 
flavour conversion. The physical argument to support the 
mean-field assumption is that the MSW resonances occur 
in the outer layers of the star, because of the large neu- 
trino mass-squared differences, and that the MSW effect 
is well accounted for in the mean-field approximation. 
However, this picture - well separating the region where 
flavour conversion occurs and the very dense one where 
the mean-free path characterizes neutrino propagation - 
has been modified by recent theoretical studies. The in- 
clusion of the neutrino-neutrino interaction in the treat- 
ment of the neutrino propagation, within the mean-field 
approximation has shown that significant flavour modifi- 
cation can occur in a region close to the neutrinosphere 
with a possible impact on nucleosynthesis [58| and 
maybe on the explosion dynamics. Ref.[57| furnishes an 
example showing that a careful treatment of the transi- 
tion region might be necessary, even from the point of 
view of the neutrino emission. Therefore, although the 
first investigations show that flavour conversion occurs 
out of the relevant region to influence the shock waves 
(see e.g. [5{J), it might still be too early to draw definite 
conclusions. 

In the context of core-collapse supernovae, a signifi- 
cant step forward beyond the mean-field approximation 
in the investigation of neutrino flavour conversion would 



require the solution of the Boltzmann equation imple- 
menting neutrino masses and mixings. Such numerical 
simulations are demanding, and appear still ahead also 
within realistic supernova simulations. The situation is 
different in the context of the early universe, since colli- 
sions are an essential ingredient bringing neutrino plasma 
to equilibration, while the mixings bring the system close 
to flavour equilibration. The neutrino history is usually 
determined by solving the Boltzmann equation for parti- 
cles with mixings, or approximate versions of it. 

Several works have addressed many-body aspects of 
the problem of neutrino flavour evolution focussing either 
on the theoretical formulation and the inherent symme- 
tries [HI , or on possible implications in schematic models 
(see e.g. (35j). In the present work we have adopted a 
novel perspective with respect to previous works and em- 
ployed the BBGKY hierarchy as a theoretical framework 
to go beyond the mean-field approximation, in a realistic 
treatment of the neutrino interactions with matter and 
among themselves. In particular we have included con- 
tributions that have not been implemented so far. The 
BBGKY hierarchy formulation offers a natural trunca- 
tion scheme with respect to the order of correlations. 
Such a formulation is equivalent to a description in terms 
of Green's function in the equal time limit. In particular, 
the BBGKY hierarchy furnishes a theoretical scheme to 
go from a many-body treatment to an effective one-body 
treatment. We have included, for the first time, novel 
vT> correlations of the pairing type. We have obtained 
coupled non-linear evolution equations for the normal 
and abnormal densities, that depend upon the normal 
and pairing mean-fields. The explicit expression for such 
mean-fields has been obtained using the usual low energy 
limit for the charged- and neutral-current interactions. 

The abnormal densities involve the expectation val- 
ues of the bilinear products of the neutrino and anti- 
neutrino operators. Such pairs correspond to neutrinos 
with the same lepton number, while different individual 
lepton numbers can appear because of the presence of 
mixings, giving rise to off-diagonal contributions of the 
abnormal density. Further conditions are imposed to the 
neutrino-antineutrino pairs, in the calculation of the ab- 
normal mean-fields, depending on the specific properties 
of the background that neutrinos are traversing. Requir- 
ing that the neutrino-antineutrino pairs have opposite 
momentum corresponds to a homogeneous medium. For 
this case we have derived a pairing mean-field with dif- 
ferent helicity conditions for the pairs. In particular, if 
one considers contributions from states of positive (neg- 
ative) helicity for the neutrinos (antincutrinos) , the pair- 
ing mean-field turns out to depend on the ratio, of the 
neutrino mass over its energy, squared (as one expects). 
Although tiny, such off-diagonal contributions might give 
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non-trivial effects. More generally one could consider the 
general case of an inhomogeneous background. The price 
to pay is that one should retain two indices in momen- 
tum in the evolution equations and all the normal and 
abnormal quantities involved. 

An hypothesis made in the present work is that the 
abnormal densities are non-zero at initial time; while 
non-zero expectation values for these bilinear operator 
products can arise by implementing the collision term in 
a Boltzmann treatment. Such a collision term is not ac- 
counted for here. However we have been discussing that, 
in principle, following the BBGKY hierarchy one has the 
complete evolution equation for the two-body correlation 
function for our system of neutrinos and anti-neutrinos 
including both the collision term (also called the Born 
terms) and the terms dependent on the two-body corre- 
lation function (also called the PP terms) that have been 
assumed to depend upon the abnormal density only in 
the present work. The three-body correlation function 
contribution is expected to give higher order corrections. 
Therefore the framework discussed here could also be em- 
ployed to investigate the relevance of the linked contri- 
butions of the two-body correlation function, that have 
made the object of the present work, for the case of cos- 
mologies] neutrinos, at the epoch of Big-bang nucleosyn- 
thesis. 

Obviously, the procedure we have been describing and 
employing for our calculations is very general. While it 
has been used to derive dynamical equations for a sys- 
tem of neutrinos propagating in a medium of ordinary 
matter and neutrinos, applying the procedure to more 
general cases is straightforward. For example, one can 
consider that matter is made up of other particles, or 
that it presents of non-standard interactions. Also, our 
results can be used for an arbitrary number of neutrino 
families and, in particular, to account for the presence of 
sterile neutrinos. 

Finally, in the present manuscript, we have been fo- 
cussing on the formal aspects implied by the two-body 
correlations of the neutrino-antineutrino type, and how 
to implement them in extended evolution equations. 
Clearly numerical calculations are required to assess the 
impact of such correlations, or, more generally, of linked 
contributions of the two-body density matrix (or correla- 
tion function) on neutrino flavour conversion in a medium 
such as a core-collapse supernova, on nucleosynthesis and 
maybe the dynamics of these explosive phenomena, or 
on cosmological neutrinos at the epoch of primordial nu- 
cleosynthesis (within a Boltzmann treatment). Further 
investigations might tell us if these contributions can en- 
gender surprising features, or novel mechanisms, in this 
fascinating domain. 
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Appendix A 



For the sake of clarity we give here an explicit formu- 
lation of Eg. (fT5|) . Writing the explicit dependence of 
each quantity on the single-particle configurations, such 
an equation reads: 



lc (ik,jl) 

= [hi(p) + h 2 (p),c 1 2]( ikij i) 
+ [(1 -Pi- P2)V(1, 2)pipa(l - Pi 2 )] M 

- [(i - Pi2)pmV{\, 2)(i - Pl - p 2 )] M 

+ [(l-pi-P2)V(l,2)ci 2 ](«ji) 

- [c 12 V(l, 2)(1 - pi- p%)](iks) 

+ ^2 Wih 3), (1 - P 13 )P1C23(1 - P^)]{ik m ,jln) (A4) 



(Al) 
(A2) 

(A3) 



^[y(2,3),(l-P 23 ) /92 ci3(l-Pi2)]( i 



km,jln) 



(A5) 



where the different contributions are 
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1 = 



3 = 



^ + ri)i r C( r fc j/) — C(ik,rl) {tl + Fl) r j 
r 

+ {t-2 + ^2)krC(ir,jl) ~ c (ik ,jr) {h + ^2Jrl 

^ j^(ifc,rs)Pl,rjP2,si(l - Pjl) 
rs 

— (1 — P r s)PlArP2MsV{rs,jl) 
va—n 

— (1 — P r s)Pl,irP2,ksV(rs,ml)Pl,nj 

Pl,rjp2,sl{ 1 ~ Pjl) 

— (1 — Prs)Pl,irP2,ksV(rs,jm)P2,nl j 

^ ; ^ ^(ik,rs) C(rs,jl) C(ik,rs)t , (rs,jl) 
rs 

^ ^ j^Pl,im^(nA:,rs)C(rs j7) C(ik,rs)V(rs,ml)Pl,nj 
m—n 

+ P2,km v (in,rs) c (rs,jl) c {ik,rs) v (rs,jm)P2,nl } 



Note the minus sign that comes from the fact that there 
are anti-neutrinos in the background. Implementing that 
the expectation value over the background : 

?^ = {bl a (p,Vb Vli (S>',ti)) (B2) 

satisfies : 

Pp^ph = (2n) 3 2E p 5 hhl 5 3 (p- ?>)p7^ (B3) 



one gets : 



(p v ) = -N(k,k') 



G, 



d 3 p 



Pp 



2V2J (2n) 3 2E p 
[u Va (k, -)7 JU (1 - J5)u„ a , (k, -)] 
[%,(P,+)7 M (1-75K(P,+)] (B4) 

with N(k, k') = (2tt) 3 8 3 (k~k'). By using the well known 
trace relations : 

u Va (k, h)~i v (l ~ j 5 )u Va ,(k,h) 



Tr 
Tr 

2 (fc T; m u s r i 



u Ua , (k, h)u Ua (k, h)~f v (l - 75) 



(B5) 



with 



rs m—n 



— (1 — P r s)piSrC(km,ls)(^ — Pik)v(rsjn) 

'- ^ ^ ^(fem,rs)(l — -Prs)P2,WC(.j SJn )(l — Pj; ) 
rs m—n 

— (1 - Prs)P2,kr C (im,]s)(^ ~ Pik)v(rs,ln) 

Appendix B: Derivation of the normal mean- field 

For completeness, we write down intermediate results in 
the calculation of the normal fields for the case of neutri- 
nos interacting with antineutrinos, while the same pro- 
cedure leads to all the expressions of the normal mean- 
fields given in Section III. Let us consider the case of 
the off-diagonal contribution of T Eq. ([38]) . but associ- 
ated with neutrino of flavour v a traversing a medium of 
anti-neutrinos of a different flavour vp (FigJ2J) : 



the expression become^! 

r„ ai „,(p„) = -N(k,k>) 



G, 



d 3 p 



2V2J {2ir) 3 2E p 
4 [k v - m v s v (k)} [p v + m p s' 7 (p)] p p 

k ■ p — m v mi, s(k) ■ s(p) 

- m v s(k) ■ p + iUyk ■ s(p) pp 13 '"" , 

(B6) 

which finally gives for relativistic neutrinos the known 
relation : 



d 3 p 



d 3 p' 



T Va , v Jp v ) = -2E k ,N(k,k') V2G F 



2n) 3 2E p J (2n) 3 2E p , 
d 3 x ^{p+k-p 



d 3 p 
J2nf 



(1 - k-p)p p 
(B7) 



Gi 
2V2. 

[u„ a (k, h a )j u (l - 75)^ (k',h'p)] 
[^(p' 5 ^)7 m (1 - 75)^ a (f, h a )] 
(bl a (p,h)b Vfl (p',h')) (Bl) 



with p p °" 



22 Note that we show terms with (effective) masses while we have 
considered relativistic neutrinos and neglected such contributions 
in the calculations. 
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Appendix C: Derivation of the abnormal mean-field 

For the sake of clarity we give here more intermediate 
steps in the calculation of the abnormal mean-fields. Let 
us first consider the homogeneous case with neutrinos 
(antincutrinos) described only by negative (positive) he- 
licity eigenstates. The off-diagonal contribution to the 
pairing potential is 



A 



d 3 k 



d 3 p 



{2i:) 3 2E k J (2tt) 3 2E p 

^= [ d 3 x e i(f+k-if-k').S 
2V2J 

[u Va {k, -)7p( 1 -l<s)u Va {k',-)] 
[^(P',+)7 AI (1-75K(P,+)] 
(al(k-)bl(p,+)) (CI) 

Implementing the homogeneity condition Eq. (|57[) and 
performing a Ficrz transformation result in 



A 



v a {k' ,-),Up{-k' ,+) 



= + (2ir) 3 6 3 (p' + k>) 



Gi 



d 3 k 



2V2J {2v) 3 2E k 



[^(-fc',+)7 Al (l-75)^ Q (fc',-)] 4' 



(C2) 



The above expression can be evaluated using the fol- 
lowing definitions. We use the chiral representation. The 
four-component spinors can be expressed as 



u(k, h) = 



v(k, h) 



<E + h\k\ X h(k) 



'E-h\k\ X h(k) 



' E - h\k\ X -h{k) 
! E + h\k\ X -h{k) l 



(C3) 



(C4) 



with the two-component helicity eigenstate spinors given 
in spherical coordinates (polar angle 9 and azimuthal an- 
gle </>) by 



X-(k) = 



cos ■ 



(C5) 



which satisfy the useful relations: 
X+ (-k) = -e^x-ik) , X-(-k) = e-^ X+ (k) . (C6) 



By using the Dirac gamma matrices 
-l\ ^ /Oct" 



7° = 



1 



7 = 



-a 



> 7 = 



1 

-1 



and the Pauli spin matrices : 
'0 -1 



1 

1 



1 



1 

-1 



(C7) 



(C8) 



one obtains for the space-component of the first spinor 
product in Eq. ([C2]) : 



u Va (k, -)7(1 - 75)^ (~k, +) 
= AExHk)<JX-(-k) 
= AE ^cos 9 cos cf> — i sin <fi, cos 9 sin < 



AE 



1*1 



iky\k\ , k z ky -\- ik x \k\ 



(C9) 
~(k x + ky 

(CIO) 



while the time-component (7^ = 70) vanishes. The sec- 
ond spinor product Vp- (— k', +)7 M (1 — 75)^^ (fc', — ) gives 
simply a complex conjugate of the above result. Substi- 
tuting these expressions into the Eq. (|C2[) one obtains 
the results Eqs. ([59]) and ([65]). 

Since for the calculation of the abnormal field Eq. ([6"5]) 
we also use the following type spinor products: 



and 



SI = u Va (k, +) 7m (1 - 75)^ (-k, +) , (Cll) 



S2 = u Va (k,-)^(l-^)v Pb (-k,-) , (C12) 



where a and b refer to given neutrino mass eigenstate, we 
give the resuls for the time- and the space-components: 



-2m a 
-2m„ 



Uu a {k,+)lo{l -75)«Pi,(-fc,+) 
E a 6 

Eb k x iky 
Ea - Ikl + kl 



(C13) 
(C14) 



u v a {k,+)l{l -75)vp b (-k,+) 

= 2m a \ — < ^e~ l ^{snv9 cos <fi, sin 9 sin <fi, cos ( 
V E b 



2m„ 



I E d fe.T 



% ky k 



kl + fc2 \k\ 



(C15) 
(C16) 
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Similarly for the time and space components of S2 one 
obtains: 



u>v a (k, -)7o(l - 75)^ {-k, -) 



-2m b \ 



(k, -)7(1 - 75)«i?i -) 

= — 2m b \ — " e +l ^(sin0 cos 0, sin 6 sin 0, cos ( 
V E b 



(C17) 
(C18) 



'2lTl b \ 



k 2 y \k\ 



(C19) 
(C20) 



In the last expressions we have employed the following 
relativistic expressions for the spinors: 



u(k,-) 
u(k,+) 
v(k,+) 
v(k,-) 



f2E 



jkX-(k)\ 
X-(k) J 

Vil X + (k)J ■ 



(C21) 
(C22) 
(C23) 
(C24) 
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